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Abstract 

Distance measurements to Type la supernovae (SNe la) at cosmological distances 
indicate that the Universe is accelerating and that a large fraction of the critical 
energy density exists in a component with negative pressure. Various hypotheses on 
the nature of this "dark energy" can be tested via their prediction for the equation 
of state of this component. If the dark energy is due to a scalar field, its equation of 
state will in general vary with time and is related to the potential of the field. We 
review the intrinsic degeneracies of luminosity distance measurements and compute 
the expected accuracies that can be obtained for the equation of state parameter 
from a realistic high statistic SNe la experiment. 
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1 Introduction 

There is now strong evidence that the Universe is fiat and that matter only 
amounts to about 1/3 of the critical energy, the remaining 2/3 exhibiting a 
large and negative pressure. There are several candidates for this dark energy 
component, which can be characterised by their "equation of state", namely 
the ratio w = p/ p. For example the genuine Cosmological Constant has p = 
—p, while topological defects give p = —p/3 for strings or p = —2p/3 for 
domain walls. But dark energy could also be due to a possibly evolving scalar 
field, in which case the equation of state may vary with time (and redshift). 

Because they are performed at varying redshifts, luminosity distance mea- 
surements of Type la supernovae can, in principle, provide estimates of this 
possible varying equation of state. The luminosity distance, however, exhibits 
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strong degeneracies, that may forbid any precise determination of this equa- 
tion of state, especially if one allows it to vary with redshift. 

We first review the expressions involved in the calculation of the luminosity 
distance, and examine the main degeneracy that shows up when analysing 

a simulated (yet realistic) high statistics SNe la experiment. We then show 
how independent knowledge of Qm will limit the effects of the degeneracy 
and compute the expected statistical and systematic uncertainties that can 
be achieved with such an experiment. 



2 Basic Equations 

We assume a fiat universe made of 2 components : non-relativistic matter, 
which contributes to the critical density, and a single extra component X 
described by its equation of state as a function of redshift: 

wx{z) = px{z) / px{z) (1) 

The luminosity distance reads: 

dL{z) = {l + z)T{z) = {l + z) ^ (2) 

where r(z) is the Robertson- Walker comoving coordinate to an object seen (via 
massless photons) at a redshift z. Friedman's equation for a two component 
flat universe reads: 

H\z)^HI{pm{z)+ px{z))/p^ (3) 

where pm is the matter density, px is the dark energy density, and po is the 
total density today. The equation of state defines the way the component x 
behaves with expansion: 

px = -3i/(l+ w;x)px (4) 
Using dz — —(1 + z)Hdt, this equation can be integrated: 

pxiz) = px(0) expS r ci^'i±^4^. (5) 
Jo 1 + z' 

One may notice that Wx = —I corresponds to a constant px- For the matter 
component, pm, w = which leads to Pm{z) = pm(0)(1 -I-;^)^. Since observa- 
tions favour a negative wx ([5]), Pm/ Px very likely increases with redshift. 
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The luminosity distance can be written: 



M^) = r (6) 

^» JnH(l + Z'f + (l- fiM)px(2)/px(0) 



Wx(^) is the (unknown) function to be determined from measurements of di 
at various redshifts. For dark energy with negative pressure (wx < Oi), the 
higher the redshift, the higher the relative matter contribution in the above 
equation. 

To see how the unknown depends on the data dii one may invert Equa- 
tion 6, first by deriving it with respect to z: 



Hl{l - Qm) Jo 1 + z' 

where r{z) = diiz) / {1 + z) , and r' = dr/dz. Deriving once again yields: 

So, wx is related to the first and second derivatives of the distance in a highly 
non-linear way. Furthermore, the terms in the sums composing the numerator 
and denominator of Equation 8 are of comparable size and opposite sign: they 
typically differ by 30 % at a redshift of 1. As a consequence, small variations 
of the derivatives of r will result in large variations of wx- This is shown the 
other way around in [1] , through striking plots. 

Several papers recently addressed the problem of reconstructing wx (or equiv- 
alently the potential of a scalar field) from luminosity distance measurements. 
Some conclude that it is hopeless given expected uncertainties of large statis- 
tics SNe la measurements [1], others provide reasonably accurate estimations 
of Wx{z)j either based on Monte-Carlo experiments [2], or even on available 
SNe la data [4]. Before attempting to clarify what causes these fundamental 
differences, we will warn the reader about a potential misconception that may 
arise from Equation 8: both its numerator and denominator are insensitive to 
Hq because r scales as Hq^. Hence w depends on and not on D,mHq, nor 
on Hq. 



3 Measurements and uncertainties 



Type la supernovae, have been used to constraint cosmological parameters 
[5,6]. Using the light curve width-peak brightness empirical relation, the in- 
trinsic dispersion of peak brightness of SNe la was found to be about 0.15 
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magnitude. A measurement with an total error of 0.2 magnitude is already 
possible (it is the best resolution obtained in [5]), and we will assume conser- 
vatively this value as a standard error. It translates into relative uncertainties 
of measured luminosity distances of 10 %, which makes the current redshift 
measurement uncertainty totally neghgible. For the proposed SNAP space 
mission [7], the accessible redshift range is limited to 1.7, due mainly to the 
overwhelming integration times required to reach the expected signal to noise. 
There is also a limitation on the low redshift side due to the large area that 
has to be covered to discover a sufficient number of low redshift supernovae. 
Minimising systematic errors dictates that the whole sample is observed using 
the same apparatus. We will therefore consider as a conservative baseline that 
2000 SNe la with redshifts in the range [0.2, 1.4] can be efficiently observed 
and see how well we can reconstruct wx{z). 

If we have a dataset of SNe la with redshifts Zi,i — 1..N, we may compute 
least squares estimates of the cosmological parameters by minimising: 

where di is the measured luminosity distance, 9 is the set of cosmological pa- 
rameters to estimate, d^^Zi, 9) is the expected value of the luminosity distance 
for the redshift Zi and the cosmological parameters 6*, p is the relative error 
(assumed to be 0.1) on the measurement. 

If one assumes that wx is constant(wx = wq), 9 is a 2 component vec- 
tor (l^MjWo). To estimate a varying wx without any prejudice, one may 
parametrise wx as a polynomial: 

d 

wx = Y.^iz' (10) 

i=0 

Choosing a polynomial of z or I+t, does not make any difference from the 
estimation point of view since there is a one-to-one mapping of the coefficients. 
Using this expression for wx. Equation 6 can be rewritten (for clarity we limit 
the polynomial to d=2): 

where W{z) = (1 + ^)"'o-"'i+"'2 exp((wi - W2)z + ^2^72)- 

For models with varying Wx-, we will concentrate on estimating the 3 param- 
eters: (f2M,Wo,Wi), which is enough to illustrate where the difficulties show 
up, for 3 different universes (QmjWcWi) = (0.3,-1.0,0.0), (0.3,-0.8,0.3) and 
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(0.4,-1.0,1.0), labelled A B and C in the figures. Model A is the standard Cos- 
mological Constant, model B is close to a quintessence model with an inverse 
power law potential modified by supergravity [3], and model C is a toy model 
with a rapidly varying w. 

The (logarithmic) derivatives shown in Figure 1 determine the information 
that every supernova adds as a function of redshift. It is clear that higher 
redshifts provide more information, but this becomes less clear if one considers 
that several lower redshift objects can be measured in the time required for 
the measurement of one high redshift object. One important point to notice is 
that for a given model, all derivatives are very similar in shape which means 
that the parameter combination that a redshift probes does not vary strongly 
with redshift. 
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Fig. 1. Logarithmic derivatives of the luminosity distance w.r.t Qm, and wi, 
for the 3 models under study (sec text), as a function of redshift. One may notice 
similar shapes of the three derivatives for every model. 






To study what infiuences the variance (and covariance) of parameter estimates, 
we will restrict ourselves to the quadratic approximation of the which 
consists in linearising as a function of parameters around the chosen model. 
This is clearly not to be done for an actual estimation, where an accurate 
mapping of the actual values (or of the likelihood) is necessary to obtain 
reliable confidence contours. Within this linear approximation, the Hessian of 
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the w.r.t the parameters reads: 



19V 



N 
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F = J2hizi,9)h^{zi,9), with h{zi,9) 



1 ddL{zi,6) 



pdi 



:i2) 



All derivatives are usually evaluated at the minimum (i.e. the parameter 
estimate). As we are studying the estimation variance as a function of pa- 
rameters and dataset, we will evaluate the derivatives at the parameter value, 
which should be the average estimate value. Since 6 is vector, h is also a vec- 
tor, and F is then a matrix, called the Fisher (or information or sometimes 
weight) matrix. Every hK^ term of the sum is a matrix of rank 1 that is the 
information on the parameters that a given supernova adds. Within this linear 
approximation. The covariance matrix of the estimates is just C = F"^. 



model 


"C" 


Va 


Eigenvectors 


A 


0.2666 -0.9853 -0.9951 
-0.9853 0.2826 0.9643 

-0.9951 0.9G13 2.117 


0.00417 
0.0781 

2.1.52 


0.9418 0.3274 0.07586 
0.3126 -0.9363 0.1598 

-0.1233 0.12G8 0.9842 


B 


1.1913 -0.9996 -0.9976 
-0.999608 1.42776 0.9953 

-0.9976 0.9953 2.290 


0.0053 
0.125 

2.948 


0.8516 0.5086 0.1263 
0.3339 -0.7123 0.6172 

-0.4039 0.4834 0.7765 


C 


0.2762 -0.9964 0.6763 
-0.9964 0.5835 -0.7339 
0.6763 -0.7339 0.26285 


0.0060 
0.174 
0.675 


0.8840 0.4538 0.1111 
-0.2333 0.222 0.9465 
-0.4048 0.8627 -0.3028 



Table 1 

Covariance matrix of the estimates for our 3 models (where diagonal elements are the 
rms errors and off-diagonal elements are the correlation coefficients). Second column 
gives the square root of eigenvalues, the third gives the corresponding eigenvectors. 
One may notice that the last eigenvector totally dominates the error budget. Es- 
pecially for model B, the errors induced by the degeneracy forbid any independent 
measurement of the parameters. 

Table 1 gives C and its eigenvalues and eigenvectors within the linear approxi- 
mation for our 3 universes, and for 2000 uniformly spread supernovae over the 
redshift range, the distance to each of them measured with a 10 % accuracy 
(0.2 mag): for model B (the worst case), the ratio of major to minor axis of the 
error ellipsoid is more than 500, which means that the parameter combination 
corresponding to the last eigenvector is almost unmeasured. 

The dataset as proposed does not allow to estimate jointly f^M, "U^o ^-nd Wi 
with a decent precision. However, one should notice that the problem is very 
different if is frozen, in which case the errors on wq and Wi go down by 
more than one order of magnitude. In [2,4], Qm is fixed, unlike in [1] where 
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ilu is estimated and marginalized over. This explains partly why these papers 
reach different conclusions. 

The reason for this large degeneracy can be found by expanding in powers of 
z the expression for H'^ under the square root in Equation 11 (which is probed 
observationally through [dr / dz)~^/'^): 

= 1 + 3(1 + nxWo)z + 3/2(2 + 5 flxWo + Ox^i + 3 flxw^z^ + o{z^) 

-no 

(13) 

where Qx — ^ — ^m- 

Consider for a while Wi = 0: at first order in z, only QxWq is determined 
(as shown by the shape of the confidence contours in the ( flu, Wq) plane 
shown in [5]), the second order fixes Wq through flxWQ, which is rather weak. 
( ^xWq, ^x^l) are much closer to observable that ( Qm, Wq), but ^xWq really 
lacks physical sense. The scheme is the same for Wi: it enters through fix'W^i at 
second and third order and through QxWqWi at third order. As a consequence 
the parametrisation {flu, ^xWo, ^xWi) is not much less degenerate than the 
one considered here. 

Obviously, the main degeneracy goes away when fixing Qu- Before studying 
how the uncertainties on Qm affect w, we will evaluate the effect of optimising 
the redshift distribution of the dataset. 



4 Optimising the redshift distribution of the dataset 

To optimise the redshift distribution of the dataset, we first study the effect 
on the largest eigenvalue of C, of adding to the total sample a small number of 
supernovae (e.g. 40) in a given redshits bin. Figure 2 shows how this quantity 
evolves as a function of the redshift in which the 40 supernovae were added. 
Although no dramatic change is seen (errors scale as the plotted values). Figure 
2 shows that a small subsample collected at much higher redshift would reduce 
the errors more efficiently than collecting more events in the original redshift 
interval. 

One may want to be more radical, i.e. really optimise the redshift distribution 
within the redshift range in order to minimise this largest eigenvalue, keeping 
constant the total number of supernovae. This does not yield the same result as 
minimising the determinant of C as done in [9] where the optimal determinant 
is obtained at the cost of an even higher largest eigenvalue. However, as for 
the minimum determinant, the optimum is reached for a redshift distribution 
consisting in 3 delta functions, that are given in Table 2. The improvement 
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Fig. 2. value of the square root of the largest eigenvalue of C as a function of the 
redshift where one adds 40 supernovae to the sample. 

with respect to a flat distribution is not large enough to consider seriously such 
an extreme option, which has very severe drawbacks. Nevertheless, a moderate 
optimisation of the redshift distribution would be worth doing since it would 
lead to a gain in resolution. 



5 Expected accuracy on wq and Wi imposing an Q.u prior 

We computed in the Bayesian approach the Wo,Wi joint confidence contours, by 
marginalizing the estimate probability distribution over flu with a Gaussian 
prior, without recourse to the linear approximation. We used a Gaussian prior 
with a — 0.05, which is 2 to 3 times better than nowadays precision from 
large scale structures, and other methods. The results are shown in Figure 2. 
Despite our conservative inputs, model A and model B can be separated at 
more than 95 % CL, whatever is the true one. Model B can be separated from 
domain walls {{wq, wi) = (—2/3, 0)) with the same accuracy. 

It is interesting to come back again to the linear approximation to study how 
estimated errors scale with the (external) Qm uncertainty, and the (internal) 
precision of the luminosity distance measurements. Adding a prior is adding 
a ((JIm — ^^mO/c'")^ to X^) where Qm' is the prior value, and Qm remains 
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model 


"C" 


r,1 rrrO rrrO 
Zi ZZ ZO 


fractions 




0.1714 


—0.9720 


—0.989 


n on 
U.zO 


n Ar\ 
0.40 


A 


—0.972 


0.1941 


0.9295 


O.oz 


n /I n 
0.40 




— u.yoy 


0.9295 


1 Q /I 1 


1 /in 


n on 
O.zO 




0.7798 


— 0.99o9 


—0.9944 


n on 

u.zu 


n Q/i 


_I_> 


— 0.99o9 


0.9o9o 


0.9oo7 


n Qi 


n /1 1 
U.4i 




-0.9911 


0.9887 


1.178 


1 1 n 






0.1828 


-0.9918 


0.3483 


0.20 


0.50 


c 


-0.9918 


0.3715 


-0.4583 


0.86 


0.36 




0.3483 


-0.4583 


0.1957 


1.40 


0.14 



Table 2 

Result of optimisation of the redshift distribution to rcdiicc the largest eigenvalue 
of C (where C is given the same way as in Table 1). The 3 optimal distributions are 
rather similar. 



the parameter to be estimated. The Fisher and covariance matrix become: 



F' = F + 



fj2 _|_ pT(Jp 



(14) 



where C — and p = (1, 0, 0) is the (unit) vector describing the prior (here 
it only concerns Qu) and a is its standard error. The a posteriori variance of 
e.g. Wi (after application of the prior) reads: 



^ftall-pfa+p'V/Cii) (15) 

where the last approximation holds for « Cu (which is already the case 
with nowadays precision on Qu)- V{wi) scales with o"^ for large a, and with 
C33 for smaller ones. When a = Cii(l — Pi3)/Pi3, V{wi) is twice its minimum 
value ((7 = 0). This happens for a = 0.03,0.08,0.44 for models A,B and C 
respectively. Doing the same exercise for wq yields a — 0.05, 0.03, 0.02. With 
a = 0.05 and our assumptions for the experimental precision, the result quality 
would benefit both from a smaller a and from a better luminosity distance 
measurement. 

For an infinitely well defined Qm, we face a second "degeneracy" demonstrated 

by the elongated shape of the error contours. The orientation of the major axis 
defines a z (zmin) ior which Var{wQ+Wiz) is minimal, — 0.2 to 0.3. Models 
with similar wx{zmin) but different Wi (within ~ 0.3 in our assumptions), 
will be indistinguishable. This is shown by Figure 1 of [1] where luminosity 
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Fig. 3. 1 and 2 cr (39% filled areas, 86% dotted lines) confidence contours for {wo,wi) 
with a Gaussian prior an^^ = 0.05 (large contours), and with an infinite precision 
(small elongated contours). For model A, the 2 a contour escapes in an unphysical 
region. 

distances for different wx having similar Wx{z — 0.3) are shown to have the 
same luminosity distance behaviour. 



6 Influence of systematic errors 

The overall scale of luminosity distances (deduced from fluxes) depends on 
LHq, where L is the absolute luminosity of the candle. So distortions that 
bias the cosmological parameters are the ones that distort the shape oi dz- At 
lowest order in z, we may model: 



where a (hopefully small) accounts for example for the drift of photometric 
calibration of supernovae across the redshift range, an unknown evolution of 
supernovae, or light absorption in the intcrgalactic medium uniformly over the 
spectral range. The SNAP mission has been designed to reduce these effects as 
much as possible, and targets a precision of 0.02 mag over the redshift range 
z < 1.7 which corresponds to \a\ < 0.006. 




(16) 
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We computed the parameter biases for our 3 models, with ct^m = 0-05 and 
a = 0.01, and found neghgible shifts for (as expected), and biases for 
(wo,Wi) of (0,-0.19), (-0.01, 0.10) and (-0.02,-0.12) for models A,B and C. 
These systematic biases are all below the statistical uncertainty, and do not 
change significantly for cr^^ = 0.1 or 0.025. 



7 Conclusions 



With conservative hypotheses for the expected performances of a future high 
statistics SNe la experiment such as the proposed SNAP mission, we find 
that the equation of state parameter of the dark energy can be measured 
provided an accurate value of f^M is used as a prior. Detailed study of Type la 
supernovae such as with the planned nearby supernova factory project [8] 
could result in further reducing the peak-brightness intrinsic dispersion which 
directly affects the precision on the equation of state parameter. In addition, 
independent gain of precision can be achieved by (moderately) optimising the 
redshift distribution of the dataset. Finally, it is interesting to note that the 
required precision on flu could come from the SNAP mission itself, using 
large scale weak lensing which is sensitive to almost independently of 
dark energy. In [10], a ground based moderately deep weak lensing survey of 
5x5 degrees is estimated to measure f^M to 10% (for Qu around 0.3), reaching 
5% for a 10x10 degrees survey. 
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